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Abstract 

By using conformal-field theory, we classify the possible irrelevant operators for 
the Ising model with nearest-neighbor interactions on the square and triangular 
lattices. We analyze the existing results for the free energy and its derivatives and 
for the correlation length, showing that they are in agreement with the conformal- 
field theory predictions. Moreover, these results imply that the nonlinear scaling 
field of the energy-momentum tensor vanishes at the critical point. Several other 
peculiar cancellations are explained in terms of a number of general conjectures. 
We show that all existing results on the square and triangular lattice are consistent 
with the assumption that only nonzero spin operators are present. 



1 Introduction 



The role of the irrelevant operators in the two-dimensional Ising model with nearest- 
neighbor interactions has been extensively discussed in the literature. The first important 
result is due to Aharony and Fisher who showed, by using the exact results for the 
free energy and the magnetization in infinite volume, that the first correction to the 
susceptibility could be explained in terms of purely analytic corrections, i.e. without 
introducing any contribution due to irrelevant operators. The conclusions of Aharony 
and Fisher were strengthened by the analysis of 0, that showed that the behavior of 
X up to O(t^) was fully compatible with the absence of irrelevant operators. [| These 
results gave rise to the idea (which has never received the status of an explicit conjecture 
as far as we know, but which has been commonly accepted in the statistical-mechanics 
community) that no contribution from irrelevant operators is present in the free energy 
of the two-dimensional Ising model with nearest-neighbor interactions. Of course, such 
a statement cannot be generically correct, since the lattice Ising model shows explicit 
violations of rotational invariance that must be due to nonrotationally invariant irrelevant 
operators. In particular, in 0], from the analysis of the mass gap, irrelevant corrections 
with renormalization-group (RG) dimension y = —2 (respectively y = —4) were clearly 
identified on the square (resp. triangular) lattice. Of course, the question remained if 
these operators did contribute to the free energy. 

The analysis of the susceptibility of has been recently extended in [|, If. In 0, 
thanks to an impressive progress in the construction and analysis of the series expansions 
for the susceptibility, it was clearly shown that at least two irrelevant operators con- 
tribute to the expansion of the susceptibility for h = near the critical point. However, 
while these results show without doubts the presence of irrelevant operators, they do not 
characterize them. In particular, the identification of these irrelevant operators with the 
corresponding quasiprimary fields of the Ising Conformal Field Theory (CFT) is still an 
open problem. In this paper we try to make some progress in this direction. 

We shall address this problem in three steps: 

1] First, we shall discuss the CFT that describes the Ising model at the critical point. 
We shall list all operators that may appear as irrelevant ones in the lattice Ising 
model. 

2] Then, we shall compare the CFT predictions with the exact results for the free energy 
and for the magnetization and with the results for the susceptibility reported in 0. 
We shall see that these results are in perfect agreement with the RG and CFT, 
but have also peculiar features that can be explained if we make some additional 
hypotheses. The existence in the nearest-neighbor Ising model of exact transforma- 
tions that map the high-temperature phase onto the low-temperature one (duality 
or inversion transformations) plays here a major role, indicating that these peculiar 
features are strictly related to the (partial) solubility of the model. 

^We should also mention that recently a similarly unexpected cancellation was found in the free energy 
on the critical isotherm T — Tc . 
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3] The conclusions reached in the analysis of the infinite-volume free energy and of its 
derivatives are further strengthened by the analysis of the mass gap (exponential 
correlation length) and of the finite-size scaling of the free energy and of its thermal 
derivatives at the critical point (we use here the results of Finally, we analyze 

the finite-size scaling of the susceptibility at the critical point, showing that the 
dependence on the boundary conditions is in perfect agreement with the conjectures 
we have made. 

Since the analysis is rather involved and the reader could be lost in the technical 
details of the forthcoming sections, we anticipate here our main findings: 

• We do not find any evidence for the presence of the leading spin-zero irrelevant 
operator predicted by CFT, the energy-momentum tensor. This result was already 
anticipated in |]10|-Q for the two-dimensional square-lattice Ising model and in 



for the one-dimensional Ising quantum chain. Also, on the triangular lattice we do 
not observe the next-to-leading spin- zero irrelevant operator that has RG dimension 

y = -6. 

• As mentioned above, we find unambiguous evidence of the presence of nonzero-spin 
irrelevant operators in the spectrum. This is not surprising, since such operators 
are those that describe the lattice breaking of the rotational symmetry. What is 
surprising is that all results can be explained in terms of the following conjecture: 

"Tie only irrelevant operators which appear in the two-dimensional nearest-neighbor 
Ising model are those due to the lattice breaking of the rotational symmetry. " 

In some sense it can be considered as a renewed version of the original idea of 
Aharony and Fisher. 

Note that this conjecture applies only to the Ising model with nearest-neighbor inter- 
actions and it is not known whether other formulations of the Ising model satisfy the same 
conjecture (probably they don't!). Moreover, one must in principle distinguish between 
different lattice types. We find that both the square-lattice and the triangular-lattice 
results are compatible with the conjecture, but it remains to be understood if it may also 
hold on other less canonical lattices, for instance for honeycomb or Kagome lattices. 

This paper is organized as follows. In Sec. ^ we describe the model, set our notations, 
and report the basic results that are needed in the following analysis. In Sec. |^ we report 
the CFT analysis of the model at criticality and classify the possible irrelevant operators. 
In Sec. ^ we discuss the infinite-volume free energy and its derivatives with respect to 
/i for /i = 0. We show that the exact results and the results of have properties that 
cannot be anticipated from CFT and RG alone. In order to explain them, we put forward 



four conjectures that are justified in Sec. 4.2 on the basis of the available results. In Sec 



O, on the basis of the conjectures we have made, we obtain some general predictions for 
the susceptibility on the triangular lattice. The extension of the results of to such a 
lattice is very important in order to understand the validity of our conjectures. In Sec. 
|] we discuss the critical behavior of the exponential correlation length. The analysis on 
the triangular lattice is particularly interesting and gives strong support to the conjecture 
we have presented above. In Sec. |^ and ^ we consider the finite-size scaling of several 
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quantities at the critical point. We show that the existence of an inversion (duahty) 
transformation and the general conjecture presented above explain some peculiar features 
of the results found in . In Sec. |] we summarize the results and discuss some open 
problems. 



2 The Ising model with nearest-neighbor interactions 

The two-dimensional Ising model is defined by the partition function 

2j = g/5E(„,m) 0""0'™+'^E„ (^n^ (2.1) 

(Ti=±l 

where the spin variables (T„ are defined on the sites n of a regular lattice and take the 
values {±1}. The model has two phases: the low-temperature one, in which the Z2 
symmetry is spontaneously broken and the high-temperature one in which the symmetry 
is restored. The two phases are separated by a critical point which is located at /5 = /5c- 
In the following we will study several observables. We define^ the free-energy density 
K], the energy per site /i), the specific heat C(/?, /i), the magnetization per site 
M(/3, /i), and the susceptibility /i): 

= lim llog(Z(/3,/^)), (2.2) 

. ,,3, 

Ci,M . M. (2.4) 

M(,.A) . W). (2.5) 

.im ^ (2.6) 
In (p. 21) is the number of sites of a finite lattice. 



2.1 The square lattice 

On the square lattice 

/?c = I log (V2 + I) = 0.4406868... (2.7) 

and we will measure the deviations from the critical temperature in terms of the variable 
r introduced in [^]: 

- sinh2/5 ) . (2. 



2 Vsinh2/? 

For (3 = /3c, r = 0, while r > (resp. r < 0) for /3 < /3c (resp. /3 > /3c). 

^ Note that our definitions differ by powers of ttie temperature and by signs from the usual thermo- 
dynamic ones. This is irrelevant for our purposes. 
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We will use the exact expressions for the free-energy density and magnetization in zero 
field given by 



F(r,0) 
M(r,0) 



ilog (2cosh2 2/3) +F^''^s(r), 



where 



fc(r) 



" de , 



1+1- 



cos^ 6 



1/2 



+ T 



(2.9) 
(2.10) 

(2.11) 
(2.12) 



In this work, the duality transformation that maps the high-temperature phase onto the 
low-temperature one plays an important role. The variable r transforms naturally under 
such transformation, i.e. r — r. It is easy to verify that 



k{- 



^smg / 



k{ry 



A;(-r)-i/8(-r)-i/8M(-r,0) = A;(r)-i/V-i/«M(r, 0). 



(2.13) 

(2.14) 
(2.15) 



By using the exact expressions for the free energy and the magnetization we define two 
functions a(r) and 6(r) that will play a major role below. They are defined by requiring 



F(r,0) 
M(r,0) 



-Aa(r)2 log |a(r) I +Ao(r), 
Bh{T)\a{T)\'l\ 



where a(r), 6(r), and Aq{t) are regular functionsQof r, a(r) r for r 
A and B are constants. Explicitly we find 



a r 



r 1 + 

V 16 1536 

A;(r)l/8 ( 1 + 21^2 . 

^ ' V 128 



+ 0(r^ 



3589 



r^ + 0(r6) 



and 



Under duality. 



A 



27r' 



a -r 



-a r 



98304 
B = 2^'\ 



(2.16) 
(2.17) 

0, 6(0) = 1, and 

(2.18) 
(2.19) 

(2.20) 

(2.21) 



Although the susceptibility in zero field has not been computed exactly, its behavior for 
= 0, r — >^ is quite well known. In the asymptotic behavior of x fo^^ = in both 
phases was obtained: 



X±{r) = C'^\T\-"'k{TY/'F^{T)^Bf{r), 



(2.22) 



We will call a function regular if it has an expansion in integer powers of r for r — s- 0. 
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where F±{t) are regular functions of r, 

oo Lv^J 

^/(^) = (2.23) 

and r is defined in (|2.8| ). Here x+('r) (x-(t)) is the susceptibihty in the high- (low-) 
temperature phase. 

By a careful numerical study, reference found two additional important properties 
of F±{t). First, F±{t) are even functions of r. There is no rigorous proof, but we note 
that a similar property is satisfied by the two-point function in the large-x limit, see Sec. 
5.1 . Moreover, the results of can be written as 



F±(r) = [a(r)r-i]-^/^ [6(^)Mr)-^/«] ^ G±(a(r)), (2.24) 



where G±{z) are even functions of z, and a(r) and 6(r) are defined in Eqs. (p.l6|) , ( p.l7|) . 
Explicitly 

G±{z) = 1 - —z^ + ( /f - / + 0(z«), (2.25) 

where f^'^ are numerical coefficients reported in 0. Note the absence of the term of order 
z"^, a result that will play a major role below. 

2.2 The triangular lattice 

On the triangular lattice 

/?c = i log 3 = 0.2746531... (2.26) 

We measure the deviations from the critical temperature in terms of the variable r defined 
by 

r^'-^, (2.27) 

where v = tanh/?. Under the inversion transformation that maps the high-temperature 
phase onto the low-temperature one. 
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it transforms simply as r — — r. It is thus the analogous of the variable ( p.8| ) introduced 
ini. 

In zero field, the free-energy density is given by 

F(r, 0) = i log(4 sinh 2/5) + F'^^'^iT), (2.29) 

where 



o/ / l^l^logL^-hr - COS01 - COS</)2 - COS(0i 02)J , (2.30) 

^ 7o Jo 



^ ^ 2 
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the magnetization by (|2.10 ), where 



k{T) 



{l-vY(l + v) 



(2.31) 



Under r — > — r, relations ( 2.13|) , (p.l4|) , and (p.l5|) hold on the triangular lattice too. 

From the expressions of the magnetization and of the free energy, we can compute the 
functions a(r) and 6(r) that are defined by (|2.16|) and ( p.l7|) . In this case we obtain 



a(r) 
h{r) 
and 



T-3 47^-5 

^ ~ 24 ^ 10368 



161r 



7 



9 



1 + 



248832 
llr2 671r^ 



113191r 

1074954240 
lOllSr^ 



31791497r8 



288 165888 15925248 275188285440 



+ 0(r 



10- 



A 



2v^7r' 



B 



(2.32) 
(2.33) 

(2.34) 



As in the square-lattice case, the functions a(r) and 6(r) satisfy the duality relations 



3 Conformal field theory analysis 
3.1 Primary and secondary fields 

The Ising model at the critical point is described by the unitary minimal CFT with 



central charge c = 1/2 |T^. Its spectrum can be divided into three conformal families 
characterized by different transformation properties under the dual and Z2 symmetries of 
the model. They are the identity, spin, and energy families and are commonly denoted as 
[/], [cr], [e]. Let us discuss their features in detail. 

• Primary fields 

Each family contains an operator which is called primary field (and gives the name 
to the entire family). Their conformal weights are /i/ = 0, /ig- = 1/16 and = 1/2 
respectively. Since the RG eigenvalues are related to the conformal weights by 
y = 2 — 2h, all primary fields are relevant. 

• Secondary fields 

All the remaining operators of the three families (which are called secondary fields) 
are generated from the primary ones by applying the generators L_j and L_i of the 
Virasoro algebra defined by 

[Ln, Lm] = {n- m)Ln+m + ^2^^^'^ ~ l)'^n+m,0 • (3.1) 

It can be shown that, by applying a generator of index k, L_k or L_k, to a field 
(where (p = I,e,a depending on the case) of conformal weight /i^, a new operator 
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of weight h = + k is obtained. In general, any combination of L_j and L_j is 
allowed. If we denote with n the sum of the indices of the generators of type L_j and 
with n the sum of those of type L_j, the conformal weight of the resulting operator 
is + n + n. The corresponding RG eigenvalue is y = 2 — 2h^ — n — n. 

Nonzero spin states 

The secondary fields may have nonzero spin, which is given by the difference n — 
n. In general, one is interested in quantities that are invariant under the lattice 
rotation group, and thus in operators that belong to its identity representation. 
Since the lattice invariance group is a finite subgroup of the rotation group, in the 
lattice discretization of a scalar operator, operators that do not have spin zero, i.e. 
transform nontrivially for general rotations, may appear. The invariance group of 
the square lattice is the finite subgroup C4 (cyclic group of order four), which has 
four representations of "discrete" spin 0, 1, 2, and 3. An observable that transforms 
as a spin-j representation under the full rotation group belongs to a representation 
of discrete spin j (mod 4) under the action of C4. Therefore, a lattice scalar operator 
is expressed as a sum of continuum operators of spin 4j, j G N. Analogously, on a 
triangular lattice the rotation group is broken to the cyclic group of order six Cg. 
In this lattice scalar operator is expressed in terms of continuum operators 

of spin 6j, j G IN. 

Null vectors 

Some of the secondary fields disappear from the spectrum due to the null-vector 



conditions (see [0). In particular, this happens for one of the two states at level 2 
in the [a] and [e] families and for the unique state at level 1 in the identity family. 
From each null state one can generate, by applying the Virasoro operators, a whole 
family of null states. Hence, at level 2 in the identity family there is only one 
surviving secondary field, which can be identified with the stress-energy tensor T 
(or T). The second null vector in the a family appears at level 3 while in the e 
family it appears at level 4. This fact will play an important role in the following. 

Secondary fields generated by L_i 

Among all secondary fields, a particular role is played by those generated by the 
L_i Virasoro generator. L_i is the generator of translations on the lattice and 
as a consequence, it has zero eigenvalue on translationally invariant observables. 
Another way to state this result is that L_i can be represented as a total derivative, 
and as such it gives zero if applied to an operator which is the integral of a suitable 
density over the lattice, i.e. a translationally invariant operator. 

Quasiprimary operators. 

A quasiprimary field IQ > is a secondary field which satisfies the equation 

Li\Q>=0. (3.2) 

This condition eliminates all the secondary fields which are generated by L„i. The 
quasiprimary operators are the only ones that may appear in translationally invari- 
ant quantities. 
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3.2 Quasiprimary states and irrelevant operators. 

It is easy to construct, by using (p.2|), all the low- lying quasiprimary states. Here is the 
list of all quasiprimary operators up to level 10. 

• In the Identity family there is one quasiprimary state at levels 2, 4, and 6 and two 
quasiprimary states at levels 8 and 10; 

• In the energy family there is one quasiprimary state at levels 4, 6, 7, 8, and 9 and 
two quasiprimary states at level 10; 

• In the [cr] family there is one quasiprimary state at levels 3, 5, 6, 7, and 8 and two 
quasiprimary states at levels 9 and 10. 

For all these states it is possible to give the exact expression in terms of the Virasoro gen- 
erators (even if it becomes increasingly cumbersome as the level increases). For instance, 
in the identity family one finds 

Qi = L_2\I>, (3.3) 
Qi = (L% - |L_4) |/ >, (3.4) 

at level 2 and 4 respectively, where we have introduced the notation to denote the 
quasiprimary state at level n in the rj family. 

Let us now construct from the listed above the irrelevant operators which could 
appear in any lattice translationally invariant quantity. We list below those that have 
RG eigenvalue \y\ < 10. We will classify them by their spin, since operators of different 
spin appear on different lattices. Spin-zero operators are relevant in all cases, spin-(4?7,) 
operators appear on the square lattice, while spin-(6n) operators play a role only on the 
triangular lattice. 

The spin-0 operators are the following: 

• Identity family: 

Q2Q2 whose weight is 4 and RG eigenvalue is —2; 
QiQi whose weight is 8 and RG eigenvalue is —6; 

• Energy family: 

QIQ4 whose weight is 9 and RG eigenvalue is —7; 

• Spin family: 

QgQ^ whose weight is 6 + | and RG eigenvalue is — (4 + |); 
QsQs whose weight is 10 + | and RG eigenvalue is — (8 + |). 

On the square lattice we should consider the spin-four operators: 

• Identity family: 

Qi + Qi whose weight is 4 and RG eigenvalue is —2; 
QiQi + QiQi whose weight is 8 and RG eigenvalue is —6; 
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• Energy family: 

Qi + Q% whose weight is 5 and RG eigenvalue is —3. 

• Spin family: Q%Qj + O7Q3 whose weight is 10 + | and RG eigenvalue is — (8 + |). 
Also the spin-eight contribute on the square lattice at this order: 

• Identity family: 

Qi + Qi whose weight is 8 and RG eigenvalue is —6; 

• Energy family: 

Q% + Q% whose weight is 9 and RG eigenvalue is —7; 

• Spin family: 

Qt + Q% whose weight is 8 + | and RG eigenvalue is —(6 + 1). 
On the triangular lattice we should consider the spin-six operators: 

• Identity family: 

Qi + Qi whose weight is 6 and RG eigenvalue is —4; 
QiQi + QiQi whose weight is 10 and RG eigenvalue is —8; 

• Energy family: 

Q& + Qq whose weight is 7 and RG eigenvalue is —5; 

• Spin family: 

Qg -I- Qe whose weight is 6 -|- | and RG eigenvalue is — (4 -|- |). 

Higher-order spins contribute operators with y < —10. For instance, in the identity 
family one should consider the spin- 12 operator (5(2 + Q12 whose weight is 12 and RG 
eigenvalue is —10. 

Among these operators, the most important ones are: QiQi that has spin zero and 
y = —2 and should be considered both for the square and the triangular lattice; Ql + Ql 
(with y = —2) and Qi + Qq (with y = —4) that are the leading operators that break 
rotational invariance on the square and on the triangular lattice respectively. These 
operators can be explicitly related to the energy-momentum tensor. The relations are: 
QiQi = Tf, Qi + Qi = T'^ + f^, + Q^Q ^ + f^. These operators will play an 
important role in the following discussion. 

As a general remark, it is important to notice that, since only even-spin operators are 
of interest, the dimensions y of the operators satisfy the following conditions: y E 2Z for 
the identity family, y e 2Z -|- 1 for the energy family, and y e 2Z — | for the spin family. 

Finally, we want to discuss the role of the symmetries. On the lattice there are two 
exact symmetries that will play an important role. 
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Z2 symmetry: {h —h). Under this transformation the operators belonging to 
the identity and to the energy family are even, while the operators belonging to the 
spin family are odd. 

duality (inversion) symmetry for h = 0. This transformation maps the high- 
temperature phase onto the low-temperature one and with our choice of variable 
r (see (|2.8|) and ( p.27|) for the square and the triangular lattice respectively) it cor- 



responds to the mapping r — r. Under this transformation (see, e.g.. Appendix 
E of |T^) the identity operators are even, the energy operators are odd, while the 



[cr]-family operators do not have a well-defined behavior. 



4 Infinite-volume zero-momentum quantities for 

h = 

In this Section, using the results of Sec. ^, we shall derive the scaling behavior of the free 
energy, magnetization, and susceptibility at h = and we will compare these results with 
the exact expressions for F{t, 0) and M(r, 0) and with the results of on the square 
lattice. We will verify that the structure of these expressions is in agreement with the 
RG predictions, although the complicated logarithmic dependence found in [|] requires 
an extension of the usual scaling expressions. Moreover, the exact results and those of ^ 
have additional properties that are specific of the lattice nearest-neighbor Ising model 
and are probably not satisfied by a generic model belonging to the Ising universality 
class. All these properties can be explained if we make some general conjectures: they 
will be presented in Sec. |0| . 

We present a general analysis for the square and the triangular lattice. In particular, 
we will show that the extension of the work of to the triangular lattice would provide 
strong support for (or rule out) our conjectures. 



4.1 Renormalization-group predictions and conjectures 



We wish now to derive the asymptotic behavior of F(r, 0), M(r, 0), and x(t, 0) by using 
the RG approach and the classification of the irrelevant operators presented in Sec. |3^ . 
We write the free energy as [0 



F{T,h) = f,{r,h) + \ut\'^y*f.^ 



\Ut 

u 



Vj/yt 



Vj/yt 



(4.1) 



where fb{T,h) is a regular function^ of r and h"^, Ut and uj are nonhnear scaling fields 
associated with the temperature and with all other operators with corresponding dimen- 
sions yt = ^ and yj. They include the nonlinear scaling field associated with the magnetic 

'^Sometimes it is assumed that the bulk free energy depends on the temperature only [p^ , [l9[| . However, 
this conjecture is inconsistent with the rigorous results available for x- See [pO| for a critical discussion. 



10 



field with dimension yh = 15/8 and those associated with all irrelevant operators. Note 
the presence of the logarithmic term due to a resonance^ between the thermal and the 
identity operator which is responsible of the log-type singularity in the specific heat [17 



The nonlinear scaling fields are analytic functions of r and h that have well-defined trans- 
formation properties under h —>■ —h. Those associated with the identity and the energy 
family are even under the transformation, while those associated with the [a] family (and 
thus Uh too) are odd. For our purposes we can expand 

UtiT,h) = ;xi(r) + ^Ai(r) + 0(/i^), (4.2) 



u 



2 

--(r,/^) = ^^^(r) + ^X,{T) + 0{h'), (4.3) 



2 



Mj(r,/i) = hvj{T) + Oih^). (A A) 

The Z2-even operators belong to the identity and the energy family and thus, for h = 0, 
they have well-defined properties under duality: 







= 




-r) 


= -/i;(r 




-r) 





(4.5) 

In general, we expect /uj(0) 7^ 0, and therefore we can normalize these scaling fields 
by requiring fij{0) = 1. On the other hand, the energy-family scaling fields — including 
that associated with the temperature — vanish for r = and thus we normalize them by 
requiring /^j(t) ~ r. The spin-family fields are normalized by requiring Vj{0) = 1. 

Let us now present our basic conjectures that will be justified in Sec. |4.2| on the 
basis of the exact expressions for the free energy and the magnetization and of the results 
of [0. Two conjectures will be presented in different forms. The analysis reported here of 
the infinite- volume quantities gives only evidence for the weaker versions (cl) and (dO). 
Evidence for (c2) will be provided in Sec. 6, and evidence for (dl)/(d2) in Sec. ^]2[ As 
we will discuss, the analysis of x on the triangular lattice should be able to discriminate 
between (dl) and (d2). 

Let us now give the list of the conjectures: 

(a) Consider a [cr]-family operator, and let Vj{T) be the corresponding nonlinear scaling 
field for h —>■ 0, cf. ( [4.4| ). Then, either Vj{T) = 0, i.e. the corresponding operator is 
decoupled, or 

k{-r)-'\{-r) = k{ry'/%j{T). (4.6) 

Such a relation should be satisfied by ^^(t) since the corresponding operator does 
not decouple. 

(b) The functions f± and f± are even functions of the nonlinear scaling fields associated 
with the energy family. 



^ Since secondary fields belonging to a given family differ by integers, we expect additional multiple 
resonances and additional terms with higher powers of log|ut| in Eq. ([4.l[). Such higher powers have 



indeed been found in the analysis of x ■ 
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(cl) The functions f± depend only on the Z-even scahng fields. 

(c2) Stronger version of the previous one: The functions f± are constant. Such a con- 
jecture was already made by Aharony and Fisher 

(do) The nonlinear scaling field of the TT operator vanishes at the critical point: U'j^f{0, 0) = 



(dl) Stronger version of (dO): The operator TT decouples, i.e. Uj'f{T,h) = for all r 



(d2) Stronger version of (dl): The only irrelevant operators that appear in the Ising 
model are the non-rotationally invariant ones. 

We remark that these conjectures (in their stronger form) are sufficient to explain the 
existing data, but are by no means necessary. For instance, consider the three conjectures 
(d). All existing square-lattice results require only (dO). Conjectures (dl) and (d2) are 
supported by the results on the triangular lattice that will be presented in Sec. |5]^ and 
^ There we will show /^tt (''") = o{t*), which provides evidence for (dl), and /i(0) = for 
the scalar operator Q\Q\ with y = —6, which is our motivation for the conjecture (d2). 
We wish also to stress that, at least in principle, some properties may hold only on a very 
specific lattice type and thus the observed properties on the triangular lattice may not 
extend to the square-lattice case. 

Let us note that in the analysis of the scaling corrections the spin of the operator 
will play an important role. As we already mentioned in Sec. p.l| , all operators of spin 
4j (respectively 6j) appear in ( [4.1[ ) on the square (resp. triangular) lattice, j G IN. 
However, because of the rotational invariance of the critical theory, nonzero spin operators 
contribute only at second order in the Taylor expansion of the infinite- volume free energy 
in powers of Uj\ut\~'^^^^* ■ 

4.2 The square lattice 

Let us now use the exact results for F(r, 0) and M(r, 0) and the results of to provide 
evidence for the conjectures we made in the previous section. 

Setting h = in ( [4.1| ) we see that all scaling fields associated with the [a] family 
disappear. Since the dimensions of the operators belonging to the energy and to the 
identity family are integers we predict 



where /o(t) and /i(t) have a regular expansion in r. The functions /o(t) and /i(t) can 
in principle depend on the phase, but from the exact solution we know that this is not 
the case. This implies 



0. 



and h. 



F(r,/i = 0)± = /o(r)+/i(r) log|r 



(4.7) 



0({x,}) 
0({x,}) 



u {{xjY''; {x, = or) = t {{x.Y'^; {x, = or) . 



(4.8) 
(4.9) 
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Using ( p.l4| ), we find that /i(t) is even in r, a property that is certainly satisfied if the 
conjecture (b) is true, i.e. (f){{xj}) is an even function of the energy-family scaling fields. 
If this is true, the energy-family scaling fields would begin to contribute to second order. 

Let us now consider the magnetization in the low-temperature phase. From ( |4.1| ) we 
obtain (r < 0) 

Mir) = J2 \f^tfv,pki{fi,fi;'^y'')+\og\fit\ if^tt'^vMf^.fir'Yn, (4.10) 

fee[cr] fceH 

where the functions pk and pk depend only on the scaling fields of the Z2-even operators, 
and the sums are over all [cr]-family operators. Now, if yk is the dimension of an operator 
belonging to the [a] family, y^ = —1/8 + 2n, where n is an integer. Therefore, we predict 



M(r) = (-r)i/«Mo(r) + (-r)i/«Afi(r) log(-r). 



(4.11) 



where Mq{t) and Mi(r) are regular functions of r. Now, the exact solution gives Mi(r) = 
0, a property that is satisfied if the conjecture (cl) is true. Setting Mi(r) = 0, we find a 
perfect agreement with the exact result. 

However, the exact result satisfies an additional property: Using ( 2.15 ), we have 

k{-r)-'/HU-r) = Mo(r)A;(r)-i/«. (4.12) 

By using the fact that yj = 2n — | (resp. yj = 2n — 1, yj = 2n) for a [a] (resp. [e], [/]) 
family operator, n G Z, it is easy to verify that such an equation is automatically satisfied 
if the conjectures (a) and (b) are true. 

Let us consider the susceptibility. By differentiating (f4.1|) and using Eqs. (|4.8| ) and 
), we obtain 



X± 



+ pAt 



h=0 



20({x,}) + 0({x,})j +p^tYl 



ViVk\Pt\ 



-Vi-Vk 



ifcef 



dx, 



-{{xj})\pt\ (Afc - ykPkKfJ't ^) + 2pt\t(pi{xj}) log I/it 



fce[/],[e] 

+Pt\og\pt\ ^ i^ik,±{{Xj})ViVk\pt\ 



-Vi-Vk 



+fJ't^og\pt\ —{{xj})\nt\y''{\k-ykPk\tPt^) 



fce[/],[e] 



(4.13) 



where all functions depend only on the irrelevant Z2-even scaling fields through Xj = 
PjP^^\ (f) and (f) are defined in Eqs. (|4.8|) , ( |4.9|) , and 'ipik,± and ipik,± are second-order 
derivatives of f± and f± with respect to the [a] -family fields. The sums over Z2-even fields 
include only the irrelevant ones — the temperature should be excluded — while the sums 
over [(j]-fields include both the magnetic and the irrelevant ones. Since yj = — 1/8 + 2n, 
n integer, for [a] operators and yj integer for Z2-even operators, this result implies the 
expansion 



X± = |r|-^/^A±(r) + |r|-^/Mog |r|B±(r) + C(r) + D{t) log |r| 



(4.14) 
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where all functions are regular and only A± and B± depend on the phase. 

If we now use the conjecture (cl) we obtain ipik,± = 0, and therefore B±{t) = in 
agreement with the results of [||. 

By comparing (^A4|) with ( ^1221) , we find Bf{r) = C{t) + D{t) log |r|, so that 5/(r) 



should be identical in both phases, in agreement with the results of 0. However, we 
predict only a single log|r|, while in Q all powers appear. This means that our scaling 
Ansatz ( |4.1| ) is not correct: There are additional resonances that give rise to a more 
complicated logarithmic structure. 
For F±{t) we find 

h{r) = ^k{T)-'/\' {^y^'^' Yl ^^^,±i{^A)v^Vk^^t''-"^~'^' . (4.15) 

ik odd 

By using the conjectures (a) and (b), we can show that F±{t) is even in r, in agreement 
with the results of [Q. Note that the functions Aj(r) instead have no specific properties 
under r — >• — r and indeed -B/(t) contains all powers of r. 

Let us now discuss in more detail the consequences of Eqs. (|2.24|) and ( p. 251) . First, 



notice that the most important irrelevant operator of the [a] family (QsQs ) has dimension 
y = —4 — 1/8. Since yt = 2 — 1/8, it gives corrections of order r^. Thus, neglecting 
corrections of this order, we need to consider only the magnetic operator (the leading one) 
among the [aj-family contributions. Second, among the Z2-even operators, the leading 
ones are TT and + T^, both with y = —2. However, + is a spin-four operator 
and thus it may contribute to rotationally invariant quantities only to second order, i.e. 
it gives corrections of order r^. Therefore, the leading correction (of order r^) can only 
be due to TT. Accordingly we write: 

= -A (1 + 0i^?/irf + 0(r^)) , (4.16) 
Ph = B{l + phifi^,fiTf + 0{T^)), (4.17) 
^±,hh = C^{1 + ^±Mi/^?/^Tf + Oir^)) . (4.18) 

Then, since Prfi''') is an even function of r, we have for the functions G±{z) defined in 
(|2l2^ ) 

G± = 1 + (V>±,Mi - 2pM + 0i)^VTf (0) + O(z^). (4.19) 



By comparing with ( p.25[ ), we see that one of the following two conditions must be satisfied: 



either {'ip±^hhi — '^Phi + 0i) = or yU^t (0) = 0- Thus, unless a miraculous cancellation 
occurs, the absence of the z'^ term implies our conjecture (dO). 

Equation ( p.25|) implies also that at least one operator contributes to order and a 
different one at order r^. Note that it is not possible that the contribution of order is 
due to the nonlinear scaling field(s) already contributing to order r^. Indeed, if this were 
the case, the contribution 0{z^) in ( ^^.25] ) would be independent of the phase as the term 
O(z^) is.0 This result is perfectly compatible with the CFT results of Sec. ^that predict: 

1. At order r^, the spin-four operator + appears; 



^ Note that this independence does not follow from the RG expressions, since the functions ip+ and 
ip- are expected to be different. 
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2. At order r^, three operators may appear: the spin-zero operators QiQi and 
and the spin- four operator Ql + Ql- 

Note that + and Ql + Q\ have y = —2 and y = —3 respectively; however, since 
they have spin four, they may contribute only at second order, and therefore at O(r^) and 
O(r^) respectively. Finally, note that (|2.25|) is also in perfect agreement with the stronger 
conjecture (d2), that only non-rotationally invariant operators are present. In this case, 
we have an operator that starts contributing at order and a second one appearing at 
order r®. 

At higher orders, the situation becomes more involved. Beside the contributions of the 
expansion of the scaling fields appearing at lower orders, at order one must consider 
the fourth power of the nonlinear scaling field associated to + T"^. There is also a 
spin-zero operator QlQl with y = —7. However, because of the conjecture (b), we expect 
this operator to contribute only to second order and therefore starting at 0(r^^). 

It is interesting to note that, if the conjecture (dO) is true, Eqs. (|2.16|) and (|2.17|) 
provide the first terms of the expansion of fiti^) and Vh{T) in powers of r. Explicitly 

Mr) = r(l-Ar2 + 0(r^)), (4.20) 

v,{t) = Hry/' (l + + 0{r')Y (4.21) 



Such expansions already appear in |2^, but assume a very simple form in the variable r. 

Finally, let us see which informations we can obtain from -B/(t). As we already noted 
our expressions are not compatible with ( p.23| ) because of the presence of higher powers 
of logr. We assume here that our parametrization of the free energy gives the correct 
expression of -B/(t) up to terms of order r^, since at this order a log^ r appears. Under 
this assumption, we can compute the first terms in the expansion of Af(r). We compare 
the terms proportional to log |r|, writing 

3 

2/i,(r)Ai (r)0({O}) = b'^'''^^' + 0{t'). (4.22) 

q=l 

Using 0(0) = — l/(27r), this gives for At(r) 

00 

Ai(r)= A:(r)l/^5^Ai,r^ (4.23) 



k=0 



where 



Xto = -0.10163764897527987657904520338506263625548489685, 
Xti = 

= -0.000912698513043685863484370258366986546254622. (4.24) 

It remains unclear why, by factoring out the term A;(r)^/'^, the linear term in Aj(r) vanishes. 
Note that the value of \t2 is correct only if the conjecture (dO) holds. 
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4.3 The triangular lattice 

It is very interesting to extend the results of to the triangular lattice. Indeed, in this 
case it is possible to make a much stronger test of the conjectures we have made. 



First, it is easy to see that the exact results |15] for the free energy and the magne- 



tization are fully compatible with the conjectures we have made. Then, let us derive the 



behavior of the susceptibility. Equation (|4.14|) is lattice independent and it implies (apart 



from the logarithmic structure) ( p.22| ). Therefore, the expansion on the triangular lattice 
should also have the form ( ^.221 ). Also, according to conjectures (a) and (b), we expect 
F[t) to be even in r, where now r is defined in ( p.27| ): some evidence will be provided in 



Sec. |5.2| . Therefore, ( |2.24| ) should hold with G±{z) even in z. 

Finally, we wish to predict which powers of z should be absent in the expansion of 
G±{z). This depends on the operators that can appear. CFT predicts the following: 

1. At order we should consider TT; 

2. At order we should consider the spin-zero operators QiQi and 

3. At order we should consider the spin-six operator + Q^; 

4. At order we should consider the spin-zero operators QqQq, QlQli and the spin- 
six operators Ql + Ql, Ql + Ql- 

As we already mentioned, spin-six operators contribute to second order in rotationally 
invariant quantities. Moreover, we have not indicated powers of lower-order operators 
and the [e]-family operator Q\Q\ that, according to conjecture (b), should contribute 
corrections of order r^^. 

From this classification, we have the following possibilities: 

1. If TT is present, the term of order z^ should be present barring miraculous cancel- 
lations. 

2. If the conjecture (dO) is true, as on the square lattice, while the conjecture (dl) is 
false so that /irf'('r) ~ '^'^ 1 then the term of order z^ should be absent and the term 
of order should be nonvanishing. 

3. If the conjecture (dl) is valid, both terms of order z^ and 2;^ should be absent; 

4. If the stronger conjecture (d2) is true, i.e. if only non-rot at ionally invariant operators 
are present, the term of order z^ is also absent. More precisely, this cancellation 
would imply /i(0) = for QiQi, f (0) = for Q3Q3, and Htt{'^) ~ o(r^). We expect 
the term of order z^ to be nonvanishing since at this order the spin-six operator 
Qi + Qi should contribute. 

The triangular lattice is therefore a better testing ground for our conjectures. Indeed, 
the conjecture (dl) requires two coefficients to vanish, a very nontrivial fact. Moreover, 
we are able to distinguish between the conjectures (dl) and (d2). 
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5 The large-distance behavior of the two-point func- 
tion 



In this Section we will study the large-distance behavior of the two-point function on the 
square lattice, reviewing in part the results of |]T2[, and on the triangular lattice. The 
square-lattice analysis will confirm the validity of the conjecture (dO), i.e. /^Tf (0) = 0. 
Much more interesting is the analysis on the triangular lattice which will show that 
f^Tfi^) = o(r'*), thus providing strong support to the conjecture (dl). We will also find 
that the subleading corrections due to the zero-spin operator with y = —6 are absent, in 
agreement with the conjecture presented in the Introduction (conjecture (d2) of Sec. [4.1| ). 



5.1 The square lattice 

Let us now consider the large-distance behavior of the two-point function for h = 0, t > 0. 



For large |a;| it has the form |21 



G{x,y;T) = Z{t) 



dki dk2 



_^ 2tt 2n A,(A;) + M,(r)2^ 



where 



A,(A;) 

Z{r) 
MJr? 



4sin — + Asin^ 



h 

2 ' ' 2 ' 

t'/^ k{Tf/^ (1 + r^) V8 = 2{k{Tf - 



r2 - 1 



(5.1) 



(5.2) 

(5.3) 
(5.4) 



From these expressions, we can compute the angle- dependent correlation length ^ {9) de- 
fined from the large-distance behavior of the two-point function along a direction forming 
an angle 9 with the side of the lattice. We obtain 



m 



v^a(r) 



1 + cos 4^ + a(r)^ 



1 An 5 
cos A9 

320 9216 



cos 8^ +0(a(r) 



48 ' ' V3072 

(5.5) 

where a(r) is defined by Eqs. ( p.l6| ), ( |2.18| ). As already observed in [Q, this expansion 
shows the presence of a correction of order due to the leading irrelevant operator 
breaking rotational invariance. However, the interesting additional feature is that this 
term is the only one, i.e. there is no correction due to the rotationally invariant subleading 
operators |12[. This result is naturally interpreted: The correction we find is due to the 



spin-four operator + and there is no contribution due the scalar operator TT. At 
order there is scalar term, but this does not require the presence of a scalar operator: 
The angle-independent contribution can be interpreted as due to the square of the spin- 
four operator T"^ + T"^. Therefore, the result ( ^.5|) supports the conjecture (dO) and is 
compatible with the stronger ones (dl) and (d2). 

In [|12| we also analyzed the on-shell renormalization constant Z{t) and found no 
terms of order r^. We thought this to be a good indication of the absence of both TT 
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and T"^ + T"^. We now believe that this conclusion was a little bit too hasty. First, ( ^.3|) 
implies 

Z{t) = V8a(r)i/^6(r)2, (5.6) 

with no corrections to all orders. Of course, we cannot take this as an indication that all 
operators are absent. Moreover, there is also a conceptual problem: Z{t) is defined from 
the behavior of the two-point function at p = —iM{T) and thus we should consider the 
momentum-dependent nonlinear scaling fields as we did in for the second-moment 
correlation length. As we shall see in the next Section, no particular simplification occurs 
in the triangular case, and we find corrections of order to the expression ( |5.6|) . Thus, 
the observed cancellation is accidental and does not have any connection with the operator 
structure of the model. 

Finally, we present an argument to make plausible the fact that the functions F±{t) 
are even in r. If the short-distance part Bf{T) were absent, such a property would follow 
from the symmetry 

(_^)-l/4 k{-T)-'/\^{~T) = T-^l' k{T)-'/\^{T). (5.7) 

The interesting observation is that this symmetry property is satisfied by the large- 
distance expression of G{x,y]T). Indeed, using the expressions reported above we im- 
mediately verify that 

(_^)-i/4 k{-rr'/'G{x, y- -t) = T-^/' k{T)-'/'G{x, y- r). (5.8) 



5.2 The triangular lattice 

We now repeat the same analysis on the triangular lattice. The large- distance behavior of 
the two-point function along a side of the lattice was computed in Such expression 
was generalized in ||23| where it was conjectured that the large-distance behavior was 
given by the propagator of a Gaussian field on a triangular lattice, in analogy with the 
square-lattice expression. Therefore, 



where 



and 



G{x,y;r) 

Mk) = 
MAt? = 



Zir) 



dki 



2n/V3 
-27r/v^ 



dko 



A,(A;) + M,(r)2^ 



, 4 , 8 A:i V3k2 

4 cos fci cos — cos . 

3 3 2 2 

g 

- (cosh - l) (cosh Ifii + 2) , 



(r) = lAiTr'/\kiTf-lf'^A{r) + ,/A{7) + l) 



1/2 



3 

log A{t) 



y/1 — V + v'^ 



(5.9) 

(5.10) 

(5.11) 

(5.12) 
(5.13) 

(5.14) 
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The conjectured form ( p.9|) was checked in the high-temperature hmit [2^], by computing 
the expansion of G{x, y; r) in powers of P to order P^^. 
Note that, under r — > — r, we have 

A{-r) = (5.15) 

and 

M,(-r)2 = M,(r)^ (5.16) 

Z(-r)(-r)-i/^A;(-r)-i/^ = Z{T)T-'^^k{T)-'^\ (5.17) 

From the large-distance behavior of the two-point function we can obtain the angle- 
dependent correlation length C,{0) taken along a direction forming an angle 6 with a side 
of the triangles. We have, in terms of the function a(r) defined in ( p.l6| ), (|2.32| ), 



m ^ 



2a(T) 



^ a(r)^ cos 66' a(r)^ cos 6^ 
^ 6480 54432 
a(r)8 a(r)S cos 6^ a^rf cos 129 
^55987200 ^ 559872 18662400 



(5.18) 



This result provides a very strong check of the conjecture (d2) presented in the introduc- 
tion. Indeed, the first correction term appears only at order a(r)^ and is proportional to 
cos 6^. It is thus unambiguously related to the spin-six operator + T^. At order ^(t)^ 
there is also a correction term, but it is again proportional to cos 66* and thus it should 
be associated to a spin-six operator. Since no new operator appears at this order, it must 
be identified with an analytic correction due to the operator + T^. At order a(r)* a 
constant term and a cos 126^ appear, but they may be due to the square of the operator 

In conclusion, this calculation provides very strong evidence for the absence of TT, 
conjecture (dl) — more precisely it proves that fi^f = o(t'^) — and also for the conjecture 
(d2). Indeed, if (dl), but not (d2), were true, the spin- zero operator Ql + Q\ would 
contribute to order r^, giving rise to an angle-independent term proportional to a(r)®. 
The absence of such term supports the validity of (d2). 

Interestingly enough, this calculation allows the computation of the first analytic term 
in the scaling field /ii(r) that is associated with + T^. Indeed, if the conjecture (d2) is 
correct, the function a(r) given in ( p.32| ) coincides with the temperature scaling field at 
/i = up to terms of order r^, i.e. iitij) = o-ij) + O(r^). Then, we write 

m = (1 + a^^t{rY^ll{T) cosGO + 0(r«)) , (5.19) 

and fix a by requiring /ii(0) = 1. Then 

^,,{r) = l-^r' + 0{T'). (5.20) 
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Considering now the function Z{t), no particular simplification occurs and a correction 
term of order a(r)^ appears. Explicitly 



(5.21) 



As we already discussed in Sec. ^TI|, the presence of the quadratic term is probably related 
to the presence of a momentum-dependent contribution to the nonlinear scaling fields. 

Finally, we note that (|5.8| ) is also satisfied on the triangular lattice, as it may be easily 
shown by using ( |5.16| ) and ( p. 17] ). Again, this gives a plausibility argument for the fact 



that the function F{t) appearing in (2.22) is even on the triangular lattice too 



6 Finite-size scaling at the critical point 

Recently, there has been much effort in understanding the behavior of the Ising model in 
a finite box or strip of size L at the critical point h = t = 0, computing the finite-size free 
energy energy E^, specific heat Cl, and inverse mass gap C,l- The results obtained 
are the following: 

• In and 0, Jl and $,l were computed on a strip of width L for several different 
lattices: It was found that these two quantities have an expansion of the form 

lVl-u = Y.jxn (6-1) 

n=0 
^ oo 

n=0 



Note that in the expansion only even powers of L appear. Moreover, on a triangular 
lattice /i = /s = and Si = s^ = 0. 

Salas and Izmailian and Hu |^ computed /l. El, Cl for a square lattice L x M 
for fixed aspect ratio p = M/L, extending the results of [^]. They found: 

L\fL-U = f:^, (6.3) 

n=0 

n=0 

C, = ilogL + v^E. + fM^l. (6.5) 

n=0 

The specific heat has also been computed for a square lattice with Brascamp-Kunz 
boundary conditions in [^. However, in this case translation invariance is lost in 
one direction and thus we cannot apply straightforwardly the results presented here. 
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In this Section, we want to explain the general features of these results. 

In finite volume the general scaling expression can be generalized by writing (see, 
r8l|T9|,B7l|28|) 



e.g 



FiT,h;L) = MT,h) + ^W{{u,Ly^}) + ^\ogLW{{u,Ly^}), 



(6.6) 



where we assume that the bulk contribution is independent of L, or, more plausibly, that 
it depends on L only through exponentially small corrections 0,|T^, and the functions W 
and W depend on all scaling fields. Equation ( |6.6|) cannot be correct in general. Indeed, 
the results of indicate the presence of powers of log |r| in the susceptibility, which 
imply the presence of powers of logL in ( |6.6|) . These corrections should be relevant only 
if we consider derivatives of the free energy with respect to h, while here we set /i = 
from the beginning. In this particular case, ( |6.6|) should be correct. 

li h = 0, the [cr]-family scaling fields do not contribute, so that (B]0) becomes 



F(r,0;L) 



/,(r,0) + -lw^({^^.(r)L^^}) + -^logL W^({/i,(r)L^^ }), 



(6.7) 



where the scaling functions depend only on the Z2-even scaling fields. By using (|4.5| ) and 
the fact that the RG eigenvalues yj are even for the identity family and odd for the energy 
family we obtain 

Wi{fi,i-r)i-L)y^}) = WiMT)Ly^}) (6.8) 

and an analogous formula for W. Therefore, apart from the bulk contribution, even 
derivatives of F with respect to r contain only even powers of L, while odd derivatives 
contain only odd powers of L. This explains the particular structure of the results obtained 
by 1^-0 since 

dF 



Er 



Cr 



2V2 



dr 



V2El + 8 



T=0 

d^F 



(6.9) 
(6.10) 



r=0 



In particular, ( |6.10| ) explains why the odd terms in the expansion of Cl are related to 
those of the energy. 

For what concerns the logarithms, only Cl shows a logarithmic dependence, and only 
at leading order in L. This may be explained if 



(6.11) 



By using the results for the specific heat at criticality and in the infinite-volume limit we 
can compute the asymptotic behavior of W{x) for a: — > and x ^ 00. For x ^ 0, the 
results for Cl imply 

^ ' " (6.12) 



W{x) ^ —x'^ + 0{x^), 

271 

while in order to obtain the correct infinite- volume limit, we should have 

1 



W{x) 



27r' 



-X 



(l + 0(x-2)). 



(6.13) 
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These two results make natural the conjecture that 



W{x) = T^x^ (6.14) 

ZTT 



for all X. There are several consequences of these results: 

• Relation ( |6.11| ) and conjecture (cl) imply conjecture (c2), i.e. that the function / 
in ( [4.1|) is a simple constant, as originally suggested by Aharony and Fisher |jl|. If 
this is the case, the function /it('r) coincides with the function a(r). 



If (|6.14 ) is correct, we predict that in the expansion of d'^^F/dr'^"' at the critical 



point there is only one logarithmic term, with a coefficient that can be computed 
from the expansion of a(r). 



Let us now use ( |6.7|) to determine the corrections to the leading behavior. We obtain 

L'h = L^MO,0) + W{{x,}), (6.15) 

+ i^^L^W.({a;,}), (6.16) 





dr 


S'F , , 





^2 



=/i=0 jfcgu 



i_ j2 Ly'^y^w,,i{x,}), 

+ + (6.17) 

i&ll] 

where we write fijir) = /ij(0) + rfiij + |T^/i2,j, the functions Wi, and Wik depend only 
on the identity-family scaling fields through xj = /ij(0)L^J = L^^ and the constant A 
is defined by ( |2.16| ). We have also used the normalization conditions fii^i = 1 for the 
energy-family fields and /ij (0) = 1 for the identity-family fields. 

Let us now discuss which corrections should be expected. The important point is that 
here, at variance with the infinite-volume case, nonzero spin operators can contribute to 
first order. Indeed, the box breaks the rotational invariance down to the lattice invariance 
and therefore the mean value of a lattice operator that is not rotationally invariant but 
has the symmetries of the lattice is nonzero. This implies that no missing term is expected 
on the square lattice, in agreement with the exact results. Indeed, the lowest operator is 
the spin-four operator + that has y = —2 and belongs to the identity family, and 
is therefore able, alone, to give rise to all observed corrections. 

On the triangular lattice instead simplifications are expected. Consider first, the free 
energy f^. The absence of the term proportional to L~^, i.e. /i = 0, implies f^Tri^) = 0, 
confirming once again the conjecture (dO). The next-to- leading operator belonging to 
the identity family is the spin-six + that has y = —4. Therefore, in ( 3.15 ) the 



rp3 _j_y3 gjygg j^jgg correctlous of order L~^". The absence of the 1/L® term requires an 
additional cancellation, i.e. /i(0) for the operator QiQi that has y = —6 and zero spin, 
thereby supporting our conjecture (d2). At order 1/L® there appears a new operator 
QiQi + QiQi that gives, together with + T^, corrections of order /,-8"~4m thus 
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indistinguishable from those of + T^. At order at least the spin-12 operator 

rp6 _j_ appga];g and therefore we expect all corrections of the form /.-lO"-^™ ]jq 
nonvanishing. 

An analogous cancellation is expected for E^. For the leading correction terms are 

jl^i,tWti{xj}) + ^fii,iWii{xj}) + ... (6.18) 

where /ii(T) is the scaling field of the spin-six operator Qq+Qq that has y = —5. Reasoning 
as before, on the basis of conjecture (dO) alone, we expect no correction of order but 
the presence of all other terms. Analogously in Cl the correction should be absent. 

The results for the correlation length show the same pattern of the free energy. The 
fact that Si = S3 = on the triangular lattice provides additional evidence for the absence 
of spin-zero operators in the theory. 

It is interesting to notice that a similar finite-size scaling analysis was performed more 
than 10 years ago for the one-dimensional Ising quantum chain which belongs to the same 
universality class of the two-dimensional Ising model (for a discussion of their connection, 
see [^). In particular, in ||30| the finite-size behavior of the free energy and of the mass 
spectrum of the model was obtained and then compared in [jl3|,|3ll] with the predictions 



of perturbed CFT (see |Q for an updated review of the subject). 

Remarkably enough, also in this case the contribution of the energy-momentum tensor 
exactly disappears and the first non-zero correction is given again by the spin-four operator 

T2 + f2 pi. 



7 Finite-size scaling of the susceptibility at t = 

In the previous section we have discussed several thermal quantities at the critical point 
and verified that the observed behavior is consistent with the RG and CFT predictions 
and the conjectures we have made. Here, we want to discuss the finite-size behavior of 
the susceptibility on the square lattice, and we will check that the correction coefficients 
depend on the shape of the domain as predicted by the spin nature of the operators. 

For this purpose we study two different finite square lattices in order to verify the 
dependence of the corrections on the domain: 

Dii^ = {(no,ni) eZ^, 0<ni,n2<M-l}, (7.1) 

D[f = {(no,ni) G < ni + ^2 < 2M- 1,0 < ni -^2 < 2M- 1} . (7.2) 

In both cases the domain is a square: the first one has boundaries that are parallel to 
the lattices axes and size L = M, while the second one is rotated by 45° and has size 
L = M\/2. We use periodic boundary conditions. For domain {A) such conditions are 
obvious, for domain {B) we identify (ni, ^2) with (ni + M, 77-2 + M) and (ni + M,n2 — M). 

7.1 Renormalization-group analysis 

The finite-size scaling behavior of the susceptibility can be derived easily, starting from 
(|6.6| ). As we already said, such an expansion misses some important corrections propor- 
tional to higher powers of log L. However, they should only be of interest if we analyzed 



23 



the asymptotic behavior of x for t — 0. Here, we consider x the critical point and 
thus such corrections should vanish. 

A simple computation gives at the critical point 



Xl(0,0) 



=^^=0 fc6[/],[e 



+J2 E Ly^'^W^kiix,}), (7.3) 

ifce [a] 

where the functions depend only on the identity-family scaling fields, Xj = fij{0)L'^^ = L^^ 
and we have used the normalization conditions Vi{0) = 1, /ij(0) = 1 for spin- and identity- 
family scaling fields. 

Since i/j = 2n — | for the [crj-family operators and yj = 2n for the identity-family 
operators, where n is an integer, we have 

i,$:L---iy.,({a:,}) = L^/^5:^, (7.4) 

jfce [cr] fc=o 

i.e. the corrections contain only even powers of L. On the square lattice we do not 
anticipate any cancellation, i.e. we expect Cfc 7^ for all k. Indeed, the leading correction 
is due to the operator + T^, which has y = —2, and thus gives rise to corrections 
involving all powers of L~^. On the triangular lattice instead we expect Ci = 0, because 
of the conjecture (dO). All other terms are expected to be nonvanishing. Indeed, the 
presence of the spin-six operator + generates terms L"^", while the presence of the 
spin-six operator Qq + Qq together with the previous one generates terms L~^~^"'. 

Let us now consider the term that contains a sum over all identity- and energy-family 
operators. We expect in this case all powers of L~^, i.e. 



1 1/7 

- Yl A,(0)L-W^,({x,}) = -5:-^. (7.5) 

fce[/],[e] A:=0 



On the square lattice we should have di = 0. Indeed, the leading energy-family scaling 
field is associated with the temperature and gives a contribution of the form 

j^XtmW,{{x,}) -U^ + ^2 + Jl + --)^ (7-6) 



and thus generates all even terms in ( [7.5|) . The odd terms in (|7.5| ) are generated by the 
identity-family operators, the leading one being + T^. It gives 

^A,(0)L-W^,({a:,}) ^ 2 {j, + + ^ + ' ' ) ^ (7-7) 

and thus generates all odd terms except the first one. Hence di = 0. Note that is 
cancellation follows from CFT alone and does not require any additional hypothesis. 
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On the triangular lattice the discussion is similar although a little more complicated. 
We predict di = d2 = = dj = 0. The condition di = Q does not require any conjecture, 
while = implies the validity of the conjecture (dO). Much more interesting is to check 
whether d^ = d^ = 0, since the vanishing of these coefficients implies Xxfi^) = 
A(0) = for the operator QlQ{. Thus, the analysis of x on the triangular lattice would 
provide some additional evidence for or rule out the conjectures (dl) and (d2). 



7.2 The transfer-matrix calculation 

From the previous discussion, we can write on the square lattice 



XL 



(0, 0) = {^o + ^ + ^)+Do + (rfo + § + + 0(L-l^/^ (7.8) 



The constant Dq is lattice and geometry independent being generated by the bulk free 
energy, and it is given by 5/(0). Explicitly: 

Do = Bf{0) ^ -0.104133245093831026452160126860473433716236727314 (7.9) 

The other constants depend on the geometry of the system and in general are expected to 
be different for the two domains (A) and (B). However, this should depend on the type of 
operator that generates them. If a term is associated with a spin- zero operator, its value 
should be identical in both geometries, while if it is the first contribution of a spin-four 
operator we expect a dependence of the form cos 46', where 6 is the angle between the 
boundaries of the domain and the lattices axes. For our specific case, since 9 = ^ we 
expect the coefficient to change sign. Therefore, we predict 

4 = c^, ct = -cf, dt = dl (7.10) 

Indeed, cq and d^ are related to the magnetic and to the thermal scaling fields that have 
both spin zero. On the other hand, ci is related to the leading identity-family operator 
with y = —2. If the conjecture (dO) is correct, this term should be due only to the spin- 
four operator + and thus, according to the previous discussion, it should differ by 
a sign in the two geometries. 

In the following we shall test the predictions (|7.10| ). For this purpose it is interesting 
to note that the constants c/q and c?^ can be predicted by using the results of P,|9|,p3|]. 
Indeed, 

XtiO)Wti{x,}) =do + ^ + 0(L-3), (7.11) 



since the leading irrelevant operator contributing to (|7.5| ) has y = —2. Now, Xt{0) is given 
in ( [4.23| ), while the leading contributions to the left-hand side can be derived from the 



energy at the critical point, since 



=h=o ^ 



^'^W,{{x,}) + 0{L~% (7.12) 
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For geometry (A), using the results of we have 

1 

L2 



Wti{x,}) = wti + ^wt2 + 0(L-^), (7.13) 



where 



wti = — 1 , ^-0.220065581798270538286514481651(7.14) 

V^^^2(O)+03(O)+^^4(O) ^ ^ 

^2(0)^3(0)^4(0)1^2(0)^ + ^3(0)^ + ^4(0)^] 

9QV2 [e2(0) + e3(0) + ^4(0)]2 

0.073073526812330794515803384757 (7.15) 



Wt2 



so that 



^ 0.022366948354353361434648349198, (7.16) 
4 ~ -0.007427021467537379563283082599. (7.17) 

Note that this calculation relies only on the RG and on the CFT classification of the 
operators, but does not make use of any of the additional conjectures. 

In order to check Eqs. ( |7.8| ) and ( |7.10| ), we performed a transfer-matrix (TM) calcu- 



lation of the susceptibility. Notice that in general it is more difficult to perform a TM 
calculation in the case in which both sizes of the lattice are finite than in the case in which 
one of them is infinite, since one has to keep into account all the eigenvalues of the TM. 



7.2.1 Numerical results 

Let us see in detail the two cases that we studied: 

• Geometry (A) 

In this case we computed x on lattices of sizes up to L = 17. In order to test our 
methods we evaluated the susceptibility in two ways, by direct differentiation of the 
free energy and by using the fluctuation-dissipation theorem, i.e. by summing the 
two-point function. The results are reported in Table |I]. By comparing the two 
columns one can estimate the size of the systematic errors. 

• Geometry (B) 

In order to study geometry (B) we used the following trick. As a flrst step, we 
performed a decimation of the lattice, i.e. every second spin was integrated out. In 
this way the number of spins is reduced by half. The price one has to pay is that the 
Hamiltonian becomes more complicated and contains, in addition to the nearest- 
neighbour interaction, a next-to-nearest neighbour and a four-point interaction. In 
the presence of an external fleld also a three-point term arises. 

However, now the axes of the decimated lattice are parallel to the axes of the torus. 
Also, the new Hamiltonian only couples neighboring time slices. Therefore, we can 
apply the same TM methods used in geometry (A). 
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Table 1: Numerical estimate of the magnetic susceptibility for geometry (A). In the second 
column we give the results obtained by differentiation of the free energy and in the third 
column those obtained by summing the time-slice two-point correlation function. 



L 


Y 
A 


Y 
A 


4 


12 181742537099 


12 18174253709876 


5 


18.092431830874 


18.09243183087397 


6 


24.959397280867 


24.95939728086672 


7 


32.740662899119 


32.74066289911872 


8 


41.402340799629 


41.40234079963127 


9 


50.915891978613 


50.91589197861391 


10 


61.256768274856 


61.25676827485805 


11 


72.403538830976 


72.40353883097585 


12 


84.337262930730 


84.33726293072681 


13 


97.041023059667 


97.04102305966430 


14 


110.49957085440 


110.4995708543933 


15 


124.69905432425 


124.6990543242478 


16 


139.62680432571 


139.6268043257091 


17 


155.27116484686 


155.2711648468523 



Our numerical results are given in Table 0. We computed the magnetic suscepti- 
bility by differentiation of the free energy. The largest lattice has M = 12, which 
corresponds to L = 16.98, and is thus completely equivalent to the largest lattice 
used in geometry (A). 

7.2.2 Analysis of the data. 

We will now use the TM data to check the theoretical predictions. We expect that the 
error induced by the error on x given in Tables |I| and ^ is small compared to that due 
to the neglected higher-order corrections in ( |7.8| ). Therefore, instead of performing a fit, 
we considered as many data points as the number of free parameters of the Ansatz, and 
then required the Ansatz to be exact for them. This gives a system of equations that is 
then solved for the free parameters. We always used consecutive values of L, i.e. Li = L, 
L2 = L — l,...,Ln = L — n + 1, where n is the number of free parameters. Errors were 
estimated from the variation of the results with the lattice size and by comparison of 
different Ansatze. 

As a preliminary test we checked that y = —2 for the leading correction to scaling. 
For this purpose we studied the Ansatz 

Xl(0, 0) = L^/^ (Co + ciL^) + Do, (7.18) 

with Co, Ci, and y as free parameters. The results are summarized in Table ^. For both 
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Table 2: Numerical result for the inverse of the magnetic susceptibility for geometry (B). 



M 


1/X 


z 


n 1 /IQfi787/l 1 c;fi7/l'?1 

u. i^tyo ( o ( ^tioo ( ^toi 


3 


0.073301790137056 


4 


0.044241139068172 


5 


0.029917172878427 


6 


0.021735601983740 


7 


0.016591966498537 


8 


0.013132015183494 


9 


0.010684547791392 


10 


0.008884576737074 


11 


0.007519096948920 


12 


0.006456674647995 



Table 3: Numerical results from the Ansatz J^ ) in geometries (A) and (B) 



Geometry (A) 


L 


Co 


Cl 


y 


12 


1.0919299 


-0.0964 


-2.102 


13 


1.0919370 


-0.0915 


-2.076 


14 


1.0919414 


-0.0881 


-2.057 


15 


1.0919441 


-0.0857 


-2.044 


16 


1.0919460 


-0.0838 


-2.034 


17 


1.0919472 


-0.0823 


-2.026 


Geometry (B) 


M 


Co 


Cl 


y 


8 


1.0919297 


0.0689 


-1.922 


9 


1.0919388 


0.0720 


-1.946 


10 


1.0919435 


0.0743 


-1.962 


11 


1.0919461 


0.0761 


-1.973 


12 


1.0919477 


0.0775 


-1.982 
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geometries the numerical result for y approaches —2 as L increases. For our largest 
lattice sizes, the deviation from —2 is about 1%. In the following analysis we shall assume 
y = -2. 

Next we analyzed the data with (|7.8|) . For geometry (A), by using the known values 
of Dq, do, and d2, we found 

= 1.09195056(4) (7.19) 
cf = -0.07914(5), (7.20) 

where the quoted uncertainties were obtained by comparing the results of the Ansatz ( [7. 8] ) 
with those obtained by adding free parameter. 

For geometry (B), by using the known value of Dq, we obtain 

Co = 1.0919504(2) (7.21) 
cf = 0.0794(4), (7.22) 
= 0.019(5). (7.23) 



Our predictions ( [7.10|) appear to be very well satisfied. Moreover, our result for cq is 
in good agreement with, although much more precise than, the estimate^ of |2^, cq = 
1.09210(11). 



If we assume = and use (|7.16| ), we obtain the more precise estimate 



= 1.0919506(2) (7.24) 
cf = 0.0790(2), (7.25) 

where the error was obtained by comparing the results with and without the parameter 
d2. 

From the above analysis we see that, within the errors, the coefficients of the 1/L^ 
correction are equal in magnitude and opposite in sign for the two geometries. Since the 
two lattices are rotated by 7r/4 this implies that this correction is completely due to the 
spin-four operator + and that the scalar operator TT is absent, in agreement with 
the conjecture (dO). 



8 Concluding remarks and open issues. 

In this paper we have discussed the presently available results for the corrections to scaling 
in the two-dimensional Ising model. We have shown that all results are in agreement with 
the RG and CFT predictions. The only missing point here is a complete analysis of the 
RG resonances and consequently an extension of the scaling forms ( [4.1D and ( |6.6D to take 
into account the logarithmic structure found in . We have also shown that the existence 
of an exact symmetry in the lattice models that maps the high-temperature phase onto the 
low-temperature one plays a very important role and explains the symmetry properties 
of the results. 

'^We report here the residt of their fit with A = 7/4, since this is the correct theoretical behavior. 



29 



However, the lattice Ising model shows also features that are not predicted by CFT 
and RG and that can be explained if some additional conjectures are made. A list of 
them is reported in Sec. [4.1| . Let us summarize the evidence we have: 

• Conjectures (a) and (b). They allow to explain the symmetry properties under 
T — > — r of the free energy and of its derivatives for h = 0. Further evidence may be 
obtained by analyzing x on the triangular lattice and checking whether the functions 
F±{t) are even in r. 

• Conjecture (cl): The functions f± do not depend on the [aj-family fields. This is 
supported by the exact known results for F(r, 0) and M(r, 0) and by the results of 
H. Further evidence is obtained from the absence of a leading logarithmic correction 



in higher-point correlation functions 10, 11 



Conjecture (c2): The functions f± are constants (this is the original conjecture 
of im). The independence of f± from the Z2-even scaling fields is supported by 
the finite-size results of |^,^ discussed in Sec. ^ The conjecture follows from this 
observation and the conjecture (cl). Conjecture (c2), together with the conjectured 
formula ( 6.14| ) can be further checked by computing the logarithmic term(s) in 



d'^F/dr" at the critical point for n > 2. 

Conjecture (dO): The nonlinear scaling field of TT vanishes at the critical point. 
On the square lattice we have ample evidence in favor of (dO), which is the only 
conjecture needed to explain the existing results. Indeed, it is supported by: 

(1) The infinite-volume results of 0. 

(2) The behavior of C,{9) discussed in Sec. [O . 

(3) The dependence of x the critical point from the boundary conditions, see 
Sec. 0. 

(4) The behavior of the two-point function at the critical point, see |p3|| . 

(5) The behavior of the free energy on the critical isotherm, see 0. 

Moreover, all triangular-lattice results are compatible with it. For these reasons, 
we believe it is more than a conjecture and it is essentially proved. It is interesting 
to notice that a similar cancellation is observed in the finite-size scaling behavior 
of the free energy and of the mass spectrum in the one-dimensional Ising quantum 
chain, see IITI 



Conjecture (dl): The operator TT is decoupled. We have evidence for the validity of 
this conjecture in the triangular-lattice Ising model. The analysis of the correlation 
length C,{0) on the triangular lattice shows that /iT"f(0) vanishes at least up to terms 
of order O(r^). There are several calculations that should be feasible and would 
add further support to the validity of (dl) on the triangular lattice: 

(1) The extension of the results of Ref. M to the triangular lattice. 
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(2) The study of the dependence on the boundary conditions of the observables 
studied in Sec. |] at the critical point on the triangular lattice. This would 
unambiguously identify the spin of the leading irrelevant operator. 

(3) The study of x the critical point on a triangular lattice. It is particularly 
important to verify whether ^3, cf. ( |7.5| ), vanishes or not. If it does, it provides 
the only available evidence for XxfiS^) = 0, and thus it would strengthen the 
conjecture. 

• Conjecture (d2): Only nonzero-spin operators are present. We have evidence for 
this conjecture on the triangular lattice. The absence of spin-zero operators beside 
TT is based on the results of Sec. |5]^ and || where we showed that the existing 
exact results require /i(0) = for the spin-zero identity-family operator QiQi with 
y = —6. The studies (1) and (2) mentioned at the previous point would further check 
the conjecture. In particular, they can verify whether v{0) = for the spin-zero 
[o"]-family operator Q^Q'^ with y = —4 — |. 

Of course, as they stand, these conjectures are just "ad hoc" prescriptions, whose only 
merit is that of providing an economical way to explain all existing results. It would be 
very important to understand if there is some symmetry argument that could explain 
them. 

There remain several open questions. First of all, one may ask whether these conjec- 
tures apply to the nearest-neighbor Ising model on any regular lattice or whether some 
of them depend on the lattice structure. Another important question is how important 
the nearest-neighbor condition is: Do some of these conjectures apply also to the Ising 
model with extended interactions? Finally, one may ask whether these cancellations are 
also observed in other models. Concerning this last question, we should mention the 
results of for the three-state Potts quantum chain, which were compared with the 



CFT predictions in [|T^]. Again, the energy- momentum tensor contribution turns out 
to be compatible with zero. However, at variance with the Ising case, there is here, at 
next-to-leading order, a clear signature of a finite-size correction due to a scalar irrelevant 
operator. Even if the Potts case is slightly different from the Ising one, since this irrelevant 
operator is actually a primary operator (more precisely is the one with conformal weight 
h = this result indicates that our conjecture (d2), if true, is specific of the Ising model 
and could be somehow related to the fact that the model is soluble on the lattice. On the 
other hand, the vanishing of the correction due to the energy-momentum tensor seems to 
be a more general phenomenon. In order to understand the validity of (dO), it would be 
interesting to extend these analyses to the generic g-state Potts model or to other specific 
values of q (for instance, to percolation). 
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